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ABSTRACT 

Path-integral for theories with degenerate vacua is investigated. The origin of 
the non Borel-summability of the perturbation theory is studied. A new prescrip- 
tion to deal with small coupling is proposed. It leads to a series, which at low orders 
and small coupling differs from the ordinary perturbative series by nonperturbative 
amount, but is Borel-summable. 
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1. Introduction 



The perturbation theory plays major roles in quantum mechanics and in quan- 
tum field theories. In the latter, it often is the only reliable calculational tool one 
has. (Cross sections in the QCD are good examples of this.) In spite of its im- 
portance, the perturbation series one obtains, in most cases, has zero convergence 
radius, due to the presence of the Dyson singularities in the complex coupling- 
constant plane. 

Several years ago, many important works appeared on the convergence property 
of the perturbation series!^' As a result, the perturbation series for many models 
were found to be asymptotic series. 

If the series is asymptotic, it can be summed by using the Borel summation 
method. An alternative the Pade approximation, which is known to give conver- 
gent results for asymptotic series!^' The existence of these summation methods 
guarantees the validity of using the first few orders of the perturbation series, as 
often done in quantum field theories, such as QED. 

The situation is completely different when the theory allows tunneling. In such 
a theory, the perturbative series is simply divergent (not Borel-summable) . It has 

a cut in the physical region of the coupling constant!^ Thus the nonperturbative 
contribution of the tunneling, which in Euclidean path-integral formalism manifests 
as instantons, needs to be accounted for. Although much progress was made in 

[7] 

this approach, it is not clear how to deal with the divergence of the perturbation 
series to obtain quantitatively reliable results. 

In this paper, we will analyse this problem under a new light. Namely, in the 
path-integral formalism with finite periodicity in Euclidean time, we will argue 
that the path-integral has a singularity at the zero coupling due to the multiplicity 
of the classical minima, which does not exist in theories without tunneling. Based 
on this, we propose to re-define the path-integral so that it is equal to the original 
path-integral for small coupling but is free of the singularity. This path-integral 



2 



is suitable for the evaluation of the original path-integral for small but non-zero 
coupling. This procedure is presented in the chapter 2. The next chapter gives the 
demonstration of the power of this method for a simple one-dimensional integral. 
The chapter 4 gives the analysis of a quantum mechanical model. (Part of this 
analysis was presented in Ref.8 by one of the authors.) The last chapter gives the 
discussion and the summary. Some numerical results, which shows the validity of 
our analysis, are given in the appendix. 

2. General formalism 
2.1. The origin of non-summability 

Let us take the partition function in the Euclidean formalism. 



where g is the coupling constant of the theory, and represents all the dynamic 
variable of the theory. [Without losing generality, we assume hereafter that is 
real.] The factor A/" is a normalization factor, which we will elaborate on later. We 
also impose a periodic boundary condition in the imaginary time r e [0, (3] . This 
way, the spectrum of the theory can be obtained from the position of the poles in 
the Laplace-transform oi Z{g) with respect to /3. 

In the ordinary perturbation theory (which generates Feynman graphs in the 
quantum field theory), we first expand the integrand in powers of g; 



n=0 

The perturbative coefficients Zn, which we define by a path- integral of Fn[(j)], 




(2.1) 



oo 




(2.2) 




(2.3) 
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gives the perturbative series for the partition function; 

oo 

z^^n{g) = Y.^ng''. (2.4) 

n=0 

At the zero couphng, g — we have a free theory, and zq is expressed in terms 
of the determinant of the free theory in the usual manner. However, the small 
coupling behaviour of Z{g) is drastically different from the above: Let us denote 
the number of the degenerate minima by m. Then the functional space of has m 
points where the integrand e~'^t'^'^l is maximized. For small enough e"'^'*'^'^' has 
m well-separated "peaks" at these points. For ^ 0, peaks are completely isolated 
from other peaks and thus each peak contributes a free-theory value. Therefore we 
find that 

lim Z{g) — mzQ — mZ{Q). (2.5) 

This is the discontinuity at the zero coupling. [The situation is illustrated in Fig.l 
for m = 2.] This discontinuity is the reason for the fact that the naive perturbation 
series (2.4) is divergent (not even Borel-summable) . 

When one considers correlation functions and other observables, the existence 
of this discontinuity is obscured by implicit adjustment of the normalization factor. 
However, the discontinuity is there and the perturbative calculation leads to a 
divergent result. 

In order to illustrate these arguments, let us take a very simple analogue of 
the theory with tunneling. [Analysis of the actual quantum theory will be given 
later.] It is defined by the following "action" without any time dependence; 

S{<t^,g)^]^<t?{l-g<t^f (2.6) 

The partition function is then a simple one-dimensional integral. Obviously, for 
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zero coupling, 

/oo 
e-h^ ^AfV2^. (2.7) 
-00 

On the other hand, 

hm Z{g) = 2AfV2^, (2.8) 

because of the contribution from the two almost-Gaussian peaks at = and 
1/g. (Of course, at fl' = 0, the outer peak is at = oo and does not contribute 
to the integral. Thus (2.7).) This is the discontinuity mentioned above. In fact, 
Fig.l is the actual numerical plot for this simple model, where we have chosen that 
J\f = l/(2\/27r). A simple calculation shows that 

Zpertig) = ^ + ^ 105/ + 6930/ + (2.9) 

while for small g, Z{g) should be expressed as 

Z{g)^l + e{g) (lime(^) = 0). (2.10) 

If one is not careful about the normalization factor M ^ the difference of the first 
terms 1 and 1/2 goes unnoticed, and one would "find" the perturbative result to 
be twice as much as the (2.9); 

Z^en{g) - 1 + 6/ + 210/ + 13860/ + (2.11) 

which satisfies the property (2.10). This is what we mentioned as implicit adjust- 
ment of the normalization factor in the above. Numerically, we find that the first 
few orders of the perturbative series (2.11) (which, of course, is wrongly normal- 
ized) give excellent results for small couplings, just as perturbative series in QCD 
does. We shall see in the next chapter why this is so, in spite of the facts that the 
whole series is non-summablc and that it has wrong normalization factor. 
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2.2. A PRESCRIPTION FOR SMALL COUPLING 

The natural question is then what should be done for the small coupling ap- 
proximation. The general prescription we propose is the following. Let us assume 
the degeneracy (and tunneling) is due to a discrete symmetry of the theory* In 
such a case, the whole functional-space of (p can be divided to subspaces, which are 
mapped into each other by the symmetry. A subspace that is not mapped onto it- 
self and covers the whole space by symmetry transformations is called fundamental 
region. [This concept is commonly used in closed-string theories.] If the functional 
space is covered by m fundamental regions, the partition function (2.1) for (7 7^ 
can be written as 

Z{g)^mM [ V(f>e-^^'^^s]^ (2.12) 
Jtz 

where TZ is the fundamental region that contains the classical vacuum at = 0. 
We stress here that the above equahty does not hold for g — 0. However, since 
we are interested in the small but non-zero coupling case, we need to establish a 
small coupling approximation of the right hand side of (2.12). That expression 
would be regular at g = 0, and does not have the right value of Z{0), since in the 
limit g ^ 0, the degenerate vacua (other than the perturbative one) moves out to 
infinity and the fundamental region TZ expands to cover the whole space. In the 
simple example (2.6), the resulting expression would satisfy the property (2.10). 

In other words, the quantity defined by the right-hand side of (2.12) is regular 
at g = 0. Thus we choose to redefine the partition function by this quantity. This 
way, the observables calculated from the partition function become free from the 
zero-coupling singularity and it should make sense to construct some kind of the 
perturbation theory. 

We construct this '"improved" perturbation theory by utilizing the expansion 

★ Again, QCD is a good example of this. For cases where this is not satisfied, there are 
natural extensions of the method we present in this paper. 
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of the integrand (2.2). Namely, the small coupling approximation we obtain is 



CO „ 

Z{g) = V kn{g) kn{g) = mN / (2.13) 

Prom the previous discussion that TZ covers the whole space at zero coupling, we 
readily notice that 

lim kn{g) = mzn- (2-14) 

Therefore, for small enough coupling, the ordinary perturbation theory with nor- 
malization adjusted (such as in (2.11)) is very close to the improved perturbation 
theory. But the question is what happens for n — > oo at fixed g. In order to 
answer this question, we shall analyse the large order behaviour of Fn[(t>] in the 
next section. 

2.3. General formalism for large-order analysis 

We can express Fri[(f\ as a contour integral in the complex g^-plane as, 

where 

~S[.p,g\=S[4>,g\ + n\ogg. (2.16) 

For large n, the saddle point approximation in the complex ^f-plane is justified. 
The saddle points are defined as the solutions of 

^^[0,^] = O. (2.17) 

Since the action is usually (at most) quadratic function of g, the above saddle-point 
equation is a quadratic equation for g, whose coefficients are functional of 0. We 
shall concentrate on this case hereafter. 
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Let us now look at the choice of the integration contours. We denote the action 
as a function of g as in the following, 



S[(t>,g]^co-cig + C2g^. (2.18) 

[The coefficients co~2 are functionals of 0.] The solutions of the saddle-point equa- 
tion (2.17) are 

g^^^{l±VD), D^l-^-^. (2.19) 

4C2 

The coefficient for the Gaussian integration is 

^ ^{-D±v^){=S'^). (2.20) 



—S[(p,g] 



dg 



9=9± 2n 



We have to choose the saddle point (s) and the direction of the contour to go 
through it so that the Gaussian integration is convergent. The result is quite 
different depending on the signature of D. Namely, the whole functional space 
of (f) is divided into two parts by the signature of the functional D. We obtain 

the different expression for -F„[0] in these regions, due to the difference in the 
distribution of saddle-points. 

For > 0, we have a pair of saddle points on the real axis. Since the functional 
C2 is positive definite (as we shall see later in the examples), S" is negative at 
g-. Thus we choose the integration contour to go through the saddle point g- 
vertically to the real axis (sec Fig. 2). As a result, at the leading order of the 
Gaussian approximation, we obtain, 

1 p-S[(f>,g-] 

h-KS'Lg- 3- 



In order to find the major contribution (peak) of Fn[0], we need to solve for the 
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following equation, 



5FM dFn[ 



6(j){T) dg- 



5g- , 5Fn[0] 



(50(r) (50(r) 



= 0. (2.22) 

9- 



At the leading order of h, however, the g- derivative in the first term is zero, simply 
because of the saddle-point condition. Thus we only need to take the second term. 
This results in the ordinary equation of the motion, 

^3Mzl = 0. (2.23) 

After obtaining a solution 0(r) for a (real) it has to be substituted in the 
definition of g^ (2.19), which then becomes a self-consistent equation for g^. Only 
when one finds a real solution for this self-consistent equation, one finds a peak of 
Fn[(j)] in the D > region of the functional space of 0. 

For < 0, the saddle points g± in (2.19) are complex. We use the following 
definition. 



(We choose 9 e [0, tt], so that g^ is always in the upper half- plane.) It is convenient 
if we rewrite ci and C2 in terms of jg^j and 9 using the following relation, 

Cl = -p-r COSfc^, C2 = 7r]—F9- (2.25) 

191 ^\9r 
This way we obtain the following expressions at the saddle points g±, 

5± = CO - n - ^e^^ie, = ^sin^(sin^±icose). (2.26) 
^ 191 

Since Re(5'^^) > at both saddle points, we choose the integration contour as is 
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illustrated in Fig. 3. This results in the following approximation for 



Fn[<t>\ 



27r 



Im 



-S[(t>,g- 



(2.27) 



In this case, the equation obtained from (2.22) is, 



Im 



0, 



(2.28) 



where a is the phase 



a = arg 



1 e-^t-^'S-l 



Slgl 9- 



Using (2.24) and (2.26), we find that 



(2.29) 



(2.30) 



It is important to note here that in general the equation (2.28) is different from the 
equation of motion of the original theory. This will bring an important consequence 
later in the analysis. It also should be noted that although the values of g at the 
saddle points are complex, the resulting "effective" equation of motion 

(2.29) is real. 
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3. The simple example 



In this chapter, we shall revisit the simple example (2.6) in the previous chapter, 
using the tools developed so far. 

The coefficients Cj defined in (2.18) are; 

CO = ^0', CI = 03^ C2 = ^0^ (3.1) 

Note that due to the symmetry S{-(f), -g) = 5(0, c/), -Fn(-</') = (-1)''F„((/)). 

First, let us look at the region |0| > 2i/n, where > 0. The solutions of the 
equation of motion (2.23) are </> = 0, l/{2g-) and 1/g-. The first solution = 
corresponds io D = —oo and therefore is inconsistent. The second solution reduces 
the self-consistent equation (2.19) to the following, 

g- ^g-(^l- ^l-lQngl^ , (3.2) 

which has the solution, g- — ±l/(4y^). This corresponds io (f) — ±2^Jn, the 
boundary of the D > region. At this point, two saddle points merges, the second 
derivative S" vanishes, and the Gaussian approximation is invalid. Therefore, this 
solution is inconsistent. The third solution (f) — 1/g- leads to 

^- - Y (l - ^1 - 4^^-) ■ (3-3) 

The only real solution is g- = 0, which corresponds to = ±oo. This simply 
corresponds to the fact that is a monotonically decreasing function for large 

enough |0|. [The corresponding action is +00 and Fn{(p) = 0.] 

In the region < 2^/n {D < 0), Eq.(2.28) becomes 

0sin(7-3|^|02sin((7-e) + 2|^|Vsin((7-2e) = 0. (3.4) 
It is convenient to scale out \g\ using the new variable </? = \g\<t>- This reduces 
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Eq.(2.19) to 

\9\ = J-^ (3-5) 
V n 

and 

cos^ = — . (3.6) 

This way, (2.28) becomes equation of motion for where only parameters are n 
and 9. Thus, one could first solve it and then determine 9 self-consistently in (3.6). 
The value of \g\ is then determined from (3.5). 

In the current example, however, it is easier to obtain the solutions by com- 
bining (3.4) and (3.6) first. Then we find that the equation reduces to 

sinecos((7 + ^) = 0. (3.7) 

Among the solutions, = corresponds to the boundary of the region, D = 0, and 
is inconsistent. Thus we find n + 2 solutions labeled by an integer m = 0, ...n + 1 
as in the following, 

+ e = (m + - ) TT. (3.8) 



2 

For m <^n, the solution is approximated by 



/37r2m+l\^/^ , ^ ^ 1 /37r2m + l\^/^ 



At these points, the function F„(0) is 



The function F„(0) is the largest (and positive) for m = 0. This and the corre- 
sponding peak for m = n -|- 2 are the major contribution to the naive perturbative 
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coefficients Zn- [The solutions for m n gives the corresponding peaks for negative 
0.] The function Fn{(f)) oscillates quickly-dumped for |0| < (f)Q. The solutions (f).^ 
for n + 1 > m > 1 are the non-leading peaks. This situation can be easily seen in 
the numerical plot of Fn{4>) in Fig.4. 

After some analysis (following the spirit of Ref.4), we find that 

for even n, 

lim TV7K- 7K = , — (3.11) 

n^oo ^(j-i-l)/24ng-n/2 ^/2^ ^ ' 

[We have confirmed this result numerically for n as large as 300.] This result 
is consistent with our assertion that the peaks at m = 0, n -|- 1 are the major 
contributions. Namely, the height of these peaks reproduces the first two leading 
terms of (3.11), 

log^^n ~ Jlogn + ( 21og2 - i ) n (= login) (3.12) 



2 ° V 2 



Let us now look at the improved perturbation theory in this model. The 
discrete symmetry we mentioned before is Z2, 

(3.13) 

9 

The fundamental region we take is 

7^ = {(/)|0<^}, (3.14) 

which always contains the classical minima = 0. The multiplicity m is 2. There- 
fore, the coefficient function for the new perturbation series (2.13) is defined by 

/•l/(25) 

kn{g) = 2J\f dct>Fn{ct>). (3.15) 

J —00 

At low orders, we could look at the following expression obtained from the 
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definitions, (2.3) and (3.15), 



1 f°° 

kn{g) = 2Zn 1= / d(f)Fn{(f)) 

V2tI Jl/(2q) 



2vr Jl/{2g) 

^2zn-^ #^e-^ (3.16) 

V27r Ji/(2g) ri\ 

~ 2Zn ;= . ,„ tC 

n!V2^ (2^)3^^-! 



This expression is valid for n <^ l/^?^- Tlius tlie difference between tlie improved 
perturbation tiieory (2.13) and tiie naive perturbation tlieory (witli artificial ad- 
justment of the normalization) (2.11) is found to be non-perturbative at low orders, 
n < l/c/2. 

Let us see what happens for n » l/fl'^- We have seen in the above that the 
function Fn[(f)] has two significant peaks at ~ ±2y^. As n increases, the peak 
at ~ 2^/n moves out of the integration region < l/(2g) for n > nc= l/{16g'^). 
On the other hand, the peak at = — 2-^/n always remains in TZ. The same is true 
for lesser peaks; the positive ones moves out of 7?. for n — > oo, while the negative 
one remain. Therefore, we find that for n — > oo (with g fixed), 

kn{g) ^ {-ir~Zn. (3.17) 

Namely, the coefficient function kn{g) is positive for even n and is negative for odd 
n. Its absolute value is a smoothly increasing function of n and at even n it is half 
as much as the naive coefficient. This series is thus Borel-summable. 

We have carried out the numerical calculation of kn{g) to confirm above anal- 
ysis. The results are illustrated in Fig.5. It is apparent that the improved pertur- 
bation theory gives result that oscillates around the exact value with increasing 
amplitude, which is a typical behaviour for asymptotic series, while the naive result 
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simply diverges at higher orders. 



4. Quantum Mechanics 
We take the following action S[4>,g\, 

S[<t>, g] = dr Q^^ + l^v{g<P)^ , ^(0) = ^<f>\l - <f>f. (4.1) 

The potential is exactly the same as the action in the previous chapter and therefore 
many of the features are quite similar. This theory has an instanton solution, 

Hr) = --^ , (4.2) 

which has the action l/(6gf^). 



4.1. Evaluation of 

In this case, ci = Ss and C2 = S4/2, where 



"0 

Sm^ I dT(t>'^. (4.3) 




In the region > 0, a nontrivial solution to the classical equation of motion 
(2.23) with effective coupling is a pair of an instanton and an anti-instanton. 
This configuration was used in the analysis by Brezin ct. all*' and is worth looking 
at. Let us place the instanton at ti and the anti-instanton at T2. We assume 

-k It is interesting to note that at the transition region n ^ ric, the original cocfRcient 
Zn ^ e~-^/^^^ , since l/<i2g-^ is the value of the action at the peak </) = l/2g, which is the 
only analogue of the instanton configuration in this model. 
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that they are well-separated, 1 -C |ti — r2|(= 5t). This configuration has action 
- 1/3^2 and 



ci = ^(5r-3 + ...), ^2 = ^(5r-^ + ...). (4.4) 

(The omitted part is the interaction term between instanton and anti-instanton, 
which decreases with 5t exponentially'''' and is irrelevant here.) Using the leading 
terms, we find that (2.19) reduces to 

The only consistent solution is g- — 0, which corresponds to -F„[0] = 0. (This 
conclusion does not change when one includes the non-leading terms, as long as 
the instanton and anti-instanton are well-separated from each other.) As in the 
case (3.3) in the previous chapter, this corresponds to the uninteresting tail of 

Fn[4>\- 

One may wonder where in our analysis is the instanton pair that appeared in 
Brezin et. al.'s analysis. It turns out that it corresponds to the "wrong" saddle 
point g.^. In fact, if one replaces g- by g.^ in the above analysis, one obtains 

instead of (4.5). For 5t ^ 1, this equation has the solution gr+ = ■sj2/?,n, which is 
the coupling found in Ref.4. 

Even through g^ is the "wrong" saddle point for our purpose, their choice may 
be perfectly suitable for obtaining the perturbative coefficient Zn- The coefficient 
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Zn is written as the "double integral" 



They choose the instanton-anti-instanton pair as above and reduce the 0- integral to 
^T-integration by integrating over the nontrivial fluctuations. They integrate over 
St along the imaginary axis, and then flnally carry out the g integration through 
g+ along the real axis. This procedure was effective in calculating the Zn, but is 
useless for obtaining Fn[4>] or kn{g)- 

For D < 0, the "effective" equation of motion (2.28) can be written as follows 
by using the scaled variable (^(r) = \g\4>{T), 

8V^aV^^0_ (4.8) 



dr"^ dip 

where the effective potential Vefj(<^) is deflned by, 

1 1 

Veffif) = - '^3<^^ + 2«^4/, (4.9) 

and 

sinfcr — 6) sm(a — 29) , , „, 

«3 = ^^^, ^4= ^ . . ■ 4.10 

sm ij sm ij 

The equations (2.24) reduce to 
and 

cos. = ^, (4.12) 
where we used the following definitions, 

Rm = ^ dTip^. (4.13) 

The equations (4.11) and (4.12) corresponds to (3.5) and (3.6) in the previous 
example. The solutions can be now identified by first solving the equation of 
motion (4.8), and then by solving the self-consistent equation (4.12) for 9. 



17 



Let us look at the region ^ <^ 1, since we found the major solution in that 
region for the simple example 

in the previous chapter. We, however, assume that 

(7 ~ ^ + Ine^. (4.14) 
4 3 

is not necessarily small. Under this condition, 

^3 ~ 1 — ^ cot Cr, ~ 1 — 2^ cot (T, (4-15) 

Thus Ks,A diverges at ^ = = (97r/8n)^/^. We shall look for the solution for 
<9 <9c. Prom (4.15), we find that the effective potential (4.9) is 

Vesiv) = ^<p\l -<p){l-{l- 29 cot a)<p). (4.16) 

Therefore, for large the solution of equation of motion is a bounce solution which 
starts from (/? ~ and bounces back 

at (/9 ^ min(l//t4, 1). Using the "energy" conservation law, we can write down 
the functions R^^a as integrals 

.min(l/K4,l) 

Rm^2 ; ^ ■ 4.17 

Since ^ <^ 1, the self-consistency condition to be solved is now 

R^{k^) ^2R^{k^). (4.18) 

We have solved this equation numerically and have found a solution ^ 1.740406, 
where Ra ~ 0.143582. The value of 9 is then obtained as 

9^9,-^^-9l-rO{9l). (4.19) 

Thus we can justify the original assumption that < ^ < ^c- 
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For this solution, we obtain the following by using the virial theorem, 



co = (^'^^^S - Mi?4 ) ^ n(2 - Scoter). (4.20) 



The action is then 



11 

S -{K4 + ism2e). (4.21) 



Using this, we find that this solution yields 



1 I n \ ' «4, 



m\= , —A ^] e~^sma. (4.22) 

We have done numerical analysis of in restricted Fourier spaces. The result 

confirms that our bounce solution obtained here is in fact the position of the peak 
of Fn[4>]. The detail is given in the appendix. 

We denote this solution by (f)+ and its partner at tt — ^ by 0-(= —0+)- We 
find that 

(4.23) 

Fn[<P-] ^ (-l)"Cn-2/3n-/2A-/2, 
where C is a positive number of 0(1) and A — e~'^*/R4 ^ 1.22194. 

Since the leading behavior of the Zn is known j"' one could compare it with the 
value of Fn[(f)] at the maxima (4.23). Since -Fn[0] is an odd functional of for odd 
n, Zn is obviously zero. For even n, the known result I . Thus the most 

leading term agrees with (4.23). 
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4.2. Small g approximation 



Wc shall now discuss the improved perturbation theory. The action (4.1) is 
symmetric under the global Z2, (j){T) ^ 1/g — 0(t). The fundamental region TZ is 
then specified by the condition 

"13 



i I drm < (4.24) 



The question is then how kn{g) would behave. At the lower order of the pertur- 
bation, kn{g) differs from 2zn by only a small amount, which is non-perturbative 
in g. This can be shown as using the following relation analogous to (3.16), 

kn{g) = 2zn - 2Af [ (4.25) 

J4»l/2g 

The dominant contribution to the second term comes from the boundary of the 
integration (p = l/(2g). This term can be estimated by minimizing the free ac- 
tion under the constraint = l/{2g). Using the Lagrange-multiplier method we 
find that the solution is the constant one, = l/{2g), which has the free ac- 
tion P/{Sg^). Therefore, the dominant contribution from the second term is of 

0((/3/^2)ng-/3/(8ff^))_ 

As n increases, the second term becomes significant. The solution we have 
found in the previous section yields 

0± oc (4.26) 

where the proportionality constant is of 0(1). Therefore, as n exceeds [j3 /2gY. 
moves outside of the integration region and only 0_ remains as a major contribution 
to kn{g)- For even n, this has an effect of reducing the value of kn{g) to Zn- For 



20 



odd n, it yields negative result, which is the analytic continuation of Zn- In other 
words, for n — > cxo, the most leading terms are 



"n 




n 



_^n/2 ^ =odd. 



n/2 



for n =even 



(4.27) 



This situation is completely the same as in the simple example studied in the 
previous chapter. Our improved perturbation theory yields a Borel summable 
series. 



In this paper, wc analyzed the path-integral for theories that allow tunneling 
between degenerate classical minima, which is related by a discrete symmetry. The 
non Borel-summability of the perturbation series is explained as the result of the 
zero-coupling singularity, where the degenerate minima move out to infinity and do 
not contribute to the path-integral. For calculations at small but non-zero coupling, 
we proposed to do the path-integral only in a fundamental region of the theory. 
By calculating the large order behaviour the n-th order functional Fn[0]) we have 
identified the configuration that dominates the functional integral for the double- 
well quantum mechanical model. We have shown that half of these dominant 
configuration moves out of the fundamental region for n — >^ oo, leaving the one 
that contributes with sign (— l)'*. This shows that unlike the naive perturbation 
theory, the series expansion we propose is Borel-summable. 

Although we have carried out the actual calculation only for a double-well 
model in this paper, the formalism we developed here is applicable for a wide class 
of models, including field theoretical ones. 

In considering the calculation of the physical observables, it is important to 
consider multi-bounce configurations. In analogy with instanton configuration, 
it is natural to expect that there are multi-bounce configuration, which are not 
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exactly the solution of the equation, but is crucial in the (3 ^ oo limit. Let us 
take configuration made of b well-isolated bounces. We assume each bounce is 
the solution of (4.8). Since this configuration has R^^\ = bR^^l, (we denote the 
values for m-bounces with superscript (6)), the self- consistent equations (4.18) are 
the same as before. Thus the same value of 9, (4.19), results. Then we find 
that = Vb\g^^^\. This means that the height of the each bounce is 1/Vb 

times before. Due to this, the value of the total action is the same as before, 
sib) ^ ^(1) Therefore, for this configuration, = Although 

these configuration has lower value of -F„[0], the large measure of its collective- 
coordinate space, P^/bl, makes them important as a whole. 

Another interesting feature is the relation with the valley methods.'^" In 
these methods and in Ref.l2 in particular, the series of configurations (called "val- 
ley") that starts from the vacuum and eventually develops to well-separated in- 
stanton and anti-instanton configuration was identified for the purpose of carrying 
out the functional integral around them. The bounce-like configuration that dom- 
inates Fn,[(f)], which we found in this paper, is somewhat similar to intermediate 
configurations in this valley. This suggests a possibility that one may split the 
functional space into two parts; the one that contains the outer region of the valley 
and thus should be treated by the valley methods and the one that contains the 
lower portion of the valley. The "improved" perturbation theory for the latter 
would be Borel-summable, since the positive bounce moves out of that region into 
the non-perturbative region for n — > oo. This may be an appropriate way to treat 
the /3 — > oo case. The detailed investigation of these is in progress. 
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APPENDIX 



In this appendix, wc present the results of the numerical calculation for the 
double-well quantum mechanical model studied in the chapter 4. 

In order to evaluate numerically, we first restrict the functional space of 

0(r) to a subspace of finite dimension. We have chosen to do so by taking a finite 
number of the coefficients of the cosine Fourier-expansion of 0(r), 



1 for^^O 

(A.1) 



[2 /27r/r\ , , , ^ ^ 



All the calculation is carried out for f3 = 20, which is large enough to contain a 

[12] 

well-separated instanton-anti-instanton pair. 

First we limit ourselves to the (a°, a^) subspace. This subspace is small enough 
to allow us to calculate exact analytical expression of Fn{ao, ai) for n as large as 45 

(by using the Mathcmatica). In Fig.6 we give the plot of Fn{ao, ai) for odd n and 
even n. It is apparent that Fn{aQ,ai) has well defined peaks. We have calculated 
the position and the height of these peaks numerically. 

On the other hand, in our n — > oo asymptotic analysis, we can determine (a°, 
a^) and ^3,4 as follows: First we substitute 

00 

(^(r) = ^aV^)(r) {a'^\g\a'), (A.2) 
/=o 

into the effective equation of motion (4.8). We obtain two equations from the 
first two modes of this equation. The other two equations we need are the self- 
consistency condition (4.18) and the relation 2k,s — K4-I-I, which comes from (4.15). 
We have found only one non-trivial real positive answer, which is — 1.16625 
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and = 0.978323. This corresponds to R4 = 0.551747 and \g\ = 0.742797/x/n. 
This leads to 

(«n,as3/mp> (^n,asymp) = i'i-57008^/n, 1.31708^/1) (A.3) 

Let us compare the exact and asymptotic values. In general we expect to see 
that the exact result and the asymptotic one agree only in the leading terms, which 
is O(v^). In order to see it clearly, we consider a pair of quantities 

(C/°, = 2V^(aO - aO_i, - al_,) (A.4) 

and similarly for the asymptotic answer. This way, if the next-leading term is a 
constant, it is subtracted out of the above expression and the exact {U^, U^) should 
converge to the asymptotic value 

Kasymp, ^n,a.,..p) = (1-57008, 1.31708) (A.5) 

obtained from (A.3). 

In Fig. 7, we plot (C/°, U}^) for every fifth n as large as 45 and the point (A.5). 
We find that the convergence is in fact excellent. It is apparent that the difference 
goes to zero in proportion to 1/n. Thus the difference between a„ and an^asymp is 
of 0(nO). 

We next discuss the height of the peak of In general we assume that 

F„(a^, a^) have the form, 

Fn(aO,a^)=Cn^ni+"^A't, (A.6) 

by an analogy with that of the asymptotic answer, (4.23). Here C, A', S and 7 are 
constants. To see the convergency exphcitly, we define as F„ divided by n"^/^ 
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and consider 



Vn = (A.7) 

i'n 



If we substitute (A. 6) into the above, we obtain the following, 



(A.8) 



As n — > oo, the right hand side diverges if 7 > and goes to zero if 7 < 0. In Fig. 8, 
we find that neither happens; 7 = 0. Then the right hand side of (A.8) becomes 



Vn^A'^^+Oi-] (A.9) 



n 



for large n. On the other hand, for large n, Vn,asymp is expressed as 



, 6 



Vn,asymp ^ + O [ ^ ] (A.IO) 



The behaviour of Vn and Vn^asymp in Fig-8 is consistent with their convergence, 
which means A — A'. 

We have also investigated the peak structure in -Fio[0] in a series of subspaces 
{ao,...,am) for m = 1 ~ 4. Namely, we obtained analytical exact expression of 
Fio(ao, ■■■,am) and calculated the position of the peak numerically. The shape of 
the function (^(r) at the peak is drawn in Fig. 9 for each case. Also drawn is the 
shape of the bounce solution obtained in the main body of this paper (which, 
of course, contains all the Fourier modes). It is apparent that as the subspace is 
expanded, the shape of 0(t) approaches that of the asymptotic answer. 

Finally, we consider the five dimensional subspace (m = 4) and vary n from one 
to ten (which is the practical upper limit of our computational power). In Fig. 10 
we give the shape of the function (^(r) at the peak of for n — 1,4, 7 and 

10. We can see clearly that the approximation by the asymptotic answer becomes 
better as n gets larger. 
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In conclusion, the numerical calculation we performed strongly supports the 
dominance of the bounce solution in Fn[4>] and the validity of analytical methods 
we developed in this paper. 
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FIGURE CAPTIONS 



1) The singularity at the zero coupUng. The actual calculation is done for the 
function Z{g) for the simple model (2.6). 

2) The choice of the contour for D > (for ci < 0). 

3) The integration contour for (yf-integration for D < 0, i.e., in case there are 
two saddle-points which are complex-conjugate of each other. 

4) The behavior of Fn[(p] for the simple model for (a) n = 7 and (b) n — 8. 
Note that these are obtained by numerical calculation of the exact form of 

not from the saddle-point approximation. 

5) The value of the partition function Z ioi g = 0.1 at order n of the perturba- 
tive series. The dash-dotted line is the exact value, the dotted line the naive 
perturbation, and the solid line the improved perturbation. 

6) Plot of Fn{aQ, ai) for (a.) n — 9 and (b) n — 10. 

7) The point {U^, U^) for every fifth n from 5 to 40 is plotted. The point marked 
by a cross is the asymptotic answer (A. 5). Good convergence for n — > oo is 
apparent. 

8) Plot of Vn and Vn,asymp defined in (A. 7). These seemingly converge to a same 
value, which is a supporting evidence for the agreement of the exact result 
and the asymptotic one, 7 = and A = A'. 

9) The shape of (p{t) at the peak of -Fio[0]. The sohd hues are for the peaks 
we found in the restricted subspaces, m = 1 ~ 4, in the ascending order of 
the height of the peaks. The dotted hue is for the bounce solution we have 
found. 

10) The shape of ^^(t) at the peak of Fn[0] in the five dimensional subspace, 
m = 4. The solid lines are for n = 1,4, 7 and 10 in the ascending order of 
the height of the peaks. The dotted hue is for the asymptotic answer. We 
observe that the solid lines approach to the dotted line as n increases. 



28 



